UNCLASSIFIED 


ao2M!1« 

if  a* 

AIMED  SERVICES  TECHNICAL  INFORMAIION  ACENCT 
ARLDKTON  HALL  STAnON 
ARLINGTON  12,  VIRGINU 


UNCLASSIFIED 


RQTICE:  When  govenment  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  In  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  8. 
Qoyexnment  thereby  Incurs  no  responsibility,  nor  any 
obligation  idiatsoeirar;  and  the  feet  that  the  Oovem- 
ment  may  have  formulated,  furnished,  or  In  any  way 
supplied  the  said  drawing,  specifications,  or  other 
data  Is  not  to  be  regarded  by  InpUcatlon  or  other¬ 
wise  as  in  any  mazmer  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rl^ts 
or  permission  to  manufacture,  use  or  sell  any 
patented  Invention  that  may  In  any  way  be  related 
thereto. 


♦E.  E.  Depar+nnen+ 

University  of  New  Mexico 
Albuquerque,  New  Mexico 

This  work  was  performed  under 
Contract  Nonr  2798(01) 


ENGINEERING  EXPERIMENT  STATION 
University  of  New  Mexico 
Albuquerque,  New  Mexico 

Technical  Report  EE-7O 

Measurements  of  Electromagnetic  Atmospheric  Noise 
in  the  5-5O  cps  Region 

by 

James  F.  White* 

February  1962 

(Supplement  to  Technical  Report  EE-69) 


« 

E . E .  Department 
University  of  New  Mexico 
Albuquerque,  New  Mexico 


This  work  was  performed  under 
Contract  Nonr  2798(01) 

Reproduction  in  whole  or  in  part  is  permitted 
for  any  purpose  of  the  United  States  Government. 


CONTENTS 

Page 

List  of  Figures 
Introduction 

Experimental  Procedure  2 

Theory  g 

Observations  j^q 

Conclusions  23 

List  of  References  26 


11 


List  of  Figures 


Figure 

1  Block  diagram  of  data-collectlng  system 

2  Block  diagram  of  data-processlng  system 

3  Variation  In  Ionosphere  height  and  conductivity  8 
for  values  of  attenuation  coefficient 

4  Changes  in  mode  frequencies  due  to  Ionosphere  9 
height  and  conductivity 


5  Average  value  of  observed  field  strength  11 

during  O9OO-IOOO  (MST)  for  period  3I  July 

to  8  September  196I 

6  Average  value  of  observed  field  strength  12 

during  1000-1100  (MST)  for  period  3I  July 

to  8  September  196I 

7  Average  value  of  observed  field  strength  13 

during  1100-1200  (MST)  for  period  3I  July 

to  8  September  I96I 

8  Average  value  of  observed  field  strength  14 

during  I20O-I3OO  (msT)  for  period  5I  July 

to  8  September  I96I 

9  Average  values  of  field  strength  during  I5 

designated  hours  (mst)  for  period  3I  July 

to  8  September  196I 

10  Average  values  of  field  strength  during  I8 

designated  hours  for  period  II-I3  September  I96I 

11  Average  difference  of  field  strength  during  19 
O9OO-IOOO  (mst)  for  period  3I  July  to 

8  September  I96I 

12  Average  difference  of  field  strength  during  20 
1000-1100  (mst)  for  period  3I  July  to 

8  September  I96I 

13  Average  difference  in  field  strength  during  21 
1100-1200  (mst)  for  period  3I  July  to 


8  September  I96I 


ill 


List  of  Figures  (Con't) 


Figure 

14  Average  difference  In  field  strength 

during  120O-I5OO  (MST)  for  period  3I  July 
to  8  September  I96I 

13  Average  difference  of  field  strength 

during  O9OO-I3OO  (mst)  for  period  3I  July 
to  8  September  I96I 


Page 

22 

24 


iv 


Introduction 


A  research  of  the  literature  undertaken  In  September,  I960 
Indicated  that  while  considerable  theoretical  and  experimental 
work  had  been  done  In  atmospheric  noise  In  the  VLF  region, 
noise  data  In  the  ELF  region  was  rather  limited.  Very  little 
effort  has  been  specifically  directed  to  securing  evidence  of 
the  earth- ionosphere  cavity  resonant  frequencies  which  has  been 
predicted  by  Schumann.  As  a  consequence  this  Investigation  was 
Initiated  to  secure  atmospheric  noise  data  In  the  3  to  30  cps 
region  of  the  spectrum  and  In  particular  to  secure  evidence  of 
the  Schumann  mode  frequencies.  While  the  experiment  reported 
In  this  paper  was  In  progress  evidence  of  the  Schumann  frequen¬ 
cies  was  reported  by  several  lnvestlgators--Balser  and  Wagner 
according  to  Raemer  ^196l|,  Fitchen  et  al.  ^196lj,  and  Maple 
[1961].  The  results  of  these  Investigators  and  those  presented 
here  are  comparable. 
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Experimental  Procedure 


In  order  to  avoid  power  line  and  other  forms  of  man-made 
Interference,  portable  equipment  was  employed  to  collect  the 
noise  data  at  a  remote  location  on  the  Isleta  Indian  Reservation 
(5^°  5^'  N,  106°  34*  W  True) .  A  block  diagram  of  the  data- 
collectlng  system  Is  shown  In  Figure  1.  The  output  of  the  4.4 
meter  vertical  antenna  Is  fed  to  the  amplifier  through  a  pre¬ 
amplifier  which  serves  as  an  Impedance  matching  device.  A  2,000 
cps  low  pass  filter  Is  Incorporated  In  the  amplifier  to  eliminate 
high  frequency  signals . 

The  output  of  the  amplifier  Is  frequency  translated  by  a 
switching  technique  to  within  the  frequency  response  region  of 
a  tape  recorder  and  recorded  on  one  track  of  a  duo-channel  mag¬ 
netic  tape.  Simultaneously  the  basic  oscillator  (1,000  cps)  which 
drives  the  switch  Is  recorded  on  the  other  track  of  the  magnetic 
tape.  This  signal  Is  used  for  synchronization  during  the  play¬ 
back  mode . 

In  the  playback  mode,  shown  In  Figure  2,  the  synchronizing 
signal  drives  the  switch  to  re-translate  the  data  to  Its  original 
position  In  the  spectrum.  Variations  In  the  tape  speed  during 
record  and  playback  cause  the  same  percentage  frequency  vari¬ 
ation  In  the  data  and  synchronizing  signal.  Thus  the  data  Is 
demodulated  during  playback  In  such  a  manner  that  both  the  dat'x 
carrier  signal  and  synchronizing  signal  are  always  the  same 
frequency  and  in  phase.  Technical  Report  EE-69  (University 
New  Mexico)  describes  this  modulation  system  In  detail. 
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Fi^ur*  1.  Block  diagram  of  data-coHacting  aystom 


Figure  2.  Block  diagram  of  data -processing  system 
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After  demodulation,  the  data  from  the  awltch  le  fed  through 
a  low-paaa  filter  to  re:..jve  all  frequencies  above  30  cpe.  Ibe 
output  of  the  filter  Is  connected  to  five  filter  olianiiels 
which  are  In  parallel.  Each  channel  contains  a  selective  ampli¬ 
fier  (vftilte  Instrument  Laboratories,  Model  212)  with  a  twln-T 
filter  (White  Instrument  Laboratories,  Type  3^2) .  The  output 
of  each  selective  amplifier  Is  rectified  and  then  integrated 
by  one  channel  of  a  five-channel  analog  computer  (Oonner,  Model 
3400)  .  The  middle  20  minutes  of  each  32-mlnute  taps  Is  processed 
twice  using  five  different  filters  for  each  run.  The  ten  fre¬ 
quencies  Involved  In  this  filtering  process  are  shewn  on  the 
grajAis  of  the  data. 

In  the  complete  recording  and  recovery  process,  the  noise 
data  Is  subjected  to  two  variations  with  known  limits.  These 
variations  are  due  to  the  frequency  response  of  the  receiver 
and  changes  In  the  receiver  gain.  The  frequency  response  Is 
within  3  db  from  3  to  2000  cps .  Variation  In  gain  is  measured 
by  a  test  circuit  Incorporated  In  the  receiver  and  sMintalned 
within  1.2  db.  Under  these  limitations,  test  tapes  at  frequen¬ 
cies  In  the  3  to  30  cps  region  provide  a  means  of  determining 
the  gain  In  the  recording  and  In  the  processing  aystfm.  The 
gain  In  each  filter  channel  Is  considered  separately  during 
data  processing.  The  combined  computer  drift  and  inter-channel 
Interference  Is  less  than  2.6^.  Each  channel  has  a  .023 
unit  bandwidth.  The  output  of  the  computer  is  put  in  terms 
of  the  mean  peak  value  at  the  rectifier  Input.  Since  the  sys¬ 
tem  gain  and  antenna  effective  height  are  known  the  rectifier 
Input  is  determined  and  plotted  as  the  mean  pnak  value  of  field 
strength . 


Theory 


The  earth- ionosphere  cavity  resonant  or  mode  frequencies 
were  determined  by  Schumann  [ 1952 1  as 


n(n+l)  =  2^  f(l  +  1)  (1) 

where  r  -  radius  of  the  earth  6560  km) , 

Q 

c  -  velocity  of  light  (5  x  10  m/sec), 
and  S  -  height  of  the  ionosphere . 

This  relation  was  developed  under  the  condition  that  the 
attenuation  coefficient,  o< ,  was  equal  to  zero,  i.e.,  perfectly 
conducting  boundaries.  With  /  =  100  km,  equation (])  becomes 


^n  ""  Vn(n+1)  (2) 


which  gives  fj^  =  10.6  cps, 

f2  =  18.5  cps, 
f^  =  25.9  cps, 
=  55.5  cps. 


(5) 

etc . 


In  considering  the  conductivity  of  the 
[i960]  showed  that  the  mode  frequencies  are 
to  the  relation 


ionosphere,  Walt 
modified  according 

(4) 


where  f^  are  the  Schumann  mode  frequencies .  The  attenuation 
coefficient, oC ,  in  equation  (4)  is  given  by 
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where  M-q  =  ^  ^  10"^  h/m, 

h  -  height  of  the  Ionosphere  In  km, 
and  (T^  -  conductivity  of  the  ionosphere  in  mhos/m. 

Since  o(  is  a  function  of  h  and  <r^,  there  are  many  combina¬ 
tions  of  h  and  <r^  that  could  give  a  particular e< .  Curves  of 

equation  (5)  are  shown  in  Figure  3.  These  curves  are  normalized 

-4 

to  ot  =  .893  which  value  occurs  when  h  =  100  km  and  =  10 
mhos/m. 

With  =  0.893  ^nd  h  =  100  km,  the  Schumann  frequencies, 
f^,  given  in  (3)  are  then  modified  as  in  (4)  to  give 


A  plot  of  f^  versus  f^^  for  various  multiples  of  oC^  (normalized 
to  =  0.893)  is  shown  in  Figure  4.  Considering  h  and  to  be 
independent  of  frequency,  at  least  at  the  lower  mode  frequencies 
^Raemer  196lj ,  it  appears  that  evidence  of  a  particular  mode 
frequency  in  experimental  data  would  also  give  sufficient 
information  to  resolve  other,  perhaps  less  evident,  mode  frequen¬ 
cies  in  the  same  data.  Thus  determined  at  a  particular 
frequency  in  a  set  of  data  should  be  the  same  at  the  other  mode 
frequencies  in  that  data. 


Figure  Variation  In  ionosphere  height  and  conductivity  for 
values  of  attenuation  coefficient 


Observations 


The  data  was  recorded  during  hourly  period  from  O9OO  to 
1300  hours  MST  (1^00-2000  GMT)  dally  for  the  five  week  period, 

July  31  to  September  8,  196I.  The  observations  were  averaged 
to  obtain  an  average  field  strength  for  each  particular  hour 
during  the  complete  five  week  period.  In  addition,  the  difference 
between  the  maximum  and  mlnlmxun  observations  at  each  frequency 
for  each  hour  were  determined  on  a  weekly  basis  and  then  aver¬ 
aged  for  the  period. 

The  average  field  strengths  for  the  period  observed  during 
O9-IO,  10-11,  11-12,  amd  12-1300  are  presented  In  Figures  5,  6, 

7,  and  8,  respectively.  A  composite  presentation  for  comparison 
purposes  Is  given  In  Figure  9.  A  few  features  of  these  figures 
are  Immediately  apparent.  The  magnitude  of  the  noise  In  the  3 
to  10  cps  region  Is  approximately  twice  that  In  any  other  part 
of  the  observed  spectrum.  There  Is  a  definite  peaking  at  13.3 
cps.  The  magnitude  of  the  noise  In  the  28  cps  region  Is  approxi¬ 
mately  half  that  In  the  3  to  10  cps  region  but  greater  than 
other  parts  of  the  spectrum. 

A  study  of  these  figures  Indicates  that  the  pattern  Is 
very  similar  from  hour  to  hour.  In  addition,  the  magnitude 
of  this  pattern  Increases  with  time  during  the  O8-I3OO  dally 
observation  period.  The  peaking  and  the  symmetry  at  13. 3  cps 
is  very  consistent.  This  would  seem  to  Indicate  that  In  this 
region  the  mode  frequency  Is  very  close  to  I3.3  cps. 
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Figure  3.  Average  value  of  observed  field  strength 
during  O9OO-IOOO  (mst)  for  period 
31  July  to  8  September  I96I 
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Figure  6.  Average  value  of  observed  field  strength 
during  1000-1100  (NST)  for  period 
31  July  to  8  September  I96I 
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Figure  7.  Average  value  of  observed  field  strength 
during  1100-1200  (mst)  for  period  5I  Jbly 
to  8  Septendser  19ol 
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Figure  9.  Average  values  of  field  strength  during 

designated  hours  (MST)  for  period  5I  July 
to  8  September  1961 


VAille  the  greater  magnitude  at  2j>  cps  Is  apparent  at  some 
hours,  It  Is  not  consistently  evident.  This  fact  Is  due 
possibly  to  the  Influence  of  the  generally  Increasing  magni¬ 
tude  In  the  2j)  to  28  cps  region. 

A  study  of  the  noise  spectra  reveals  that  the  magnitude 
of  the  noise  rises  quite  sharply  In  the  lj>  to  10  cps  region. 

The  consistency  of  this  shape  seems  to  Indicate  the  existence  of 
a  peaking  In  the  7.3  to  10  cps  region.  Filter  spacing  used  In 
the  analysis  was  obviously  too  coarse  for  good  resolution  In 
this  area.  In  order  to  determine  the  existence  of  a  mode  fre¬ 
quency  In  this  region,  additional  filter  processing  could  be 
performed.  However,  In  view  of  the  pronounced  peak  at  I3.3 
cps  and  the  consistency  of  ^  for  a  set  of  data  (shown  previ¬ 
ously),  It  appears  reasonable  to  resolve  ^he  7.3  to  10  cps 
region  with  the  aid  of  equation  (6) . 

Considering  the  peak  at  I3.3  cps,  equation  (6)  becomes 


and  0^  =4.64 

or  3.2  (normalized  to  ' »  =  • 

Now  applying  3.2  o(^  to  equation  (6)  for  10.6  cps  (f^^)  gives 
f^  =  10.6  -  (3.2)(0.12o)Vl0.6 
f^  =  10.6  -  2.132  =  8.47  cps, 

23.9  cps  (f_)  gives 


for 
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and  for  55.5  cps  (fj^)  gives 

=  53.5  -  5.7856  =  29.7  cps. 

I  I 

Thus  agrees  quite  well  with  the  data.  The  frequency, 
accounts  for  the  ever-present  peaking  as  the  frequency 
Increases  to  28  cps . 

The  efficacy  of  the  above  method  In  determining  the  mode 
frequencies  can  be  demonstrated 'using  the  experimental  data 
obtained  In  this  region  of  the  spectrum  by  Balser  and  Wagner 
^1961].  Their  experimental  results  Indicate  that  f^  =  8.2  cps. 
Using  equation  (6),  e<  =  5*25  or  5.86o<q.  The  results  of  con- 

I  I  I 

sidering  In  turn  ±2  -  14.1  cps,  f^  =  20.5  cps,  and  f^^  =  26.5 
cps  show  that  they  all  involve  5.86o(q. 

The  results  of  observations  during  the  evening  hours 
(September  11  to  I5,  I96I)  are  presented  in  Figure  10.  The 
Intent  here  was  to  determine  If  variations  In  Ionosphere  height 
or  conductivity  would  be  manifested  by  a  shift  In  the  mode  fre¬ 
quencies  .  For  some  hours  there  is  a  peaking  at  18  cps .  Perusal 
of  the  data  Indicates,  however,  that  the  averages  are  peaked 
at  15.5  cps  and  thart  no  shift  is  evident.  The  25-2400  hour 
readings  which  give  the  appearance  of  I8  cps  peaking  are  some¬ 
what  erratic  and  questionable.  The  limited  amount  of  data 
precludes  noting  any  significant  shift  of  mode  frequencies 
during  the  night  hours . 

The  average  of  the  difference  between  maximum  and  minimum 
observed  field  strengths  for  the  recording  period  at  the  various 
hours  09-1500  are  shown  in  Figures  11,  12,  I5,  and  14.  The 
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Figure  11.  Average  difference  of  field  strength  during 
0900-1000  (MST)  for  period  5I  July  to 
8  September  I96I 
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Figure  12.  Average  difference  of  field  strength  during 
1000-1100  (MST)  for  period  5I  July  to 
8  September  1961 
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Figure  I5. 


Average  difference  In  field  strength 
during  1100-1200  (MST)  for  period 
51  July  to  8  September  I96I 
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Figure  14.  Average  difference  in  field  strength 
during  1200-1500  (MST)  for  period 
51  July  tc  8  September  1961 


difference  between  maximum  and  minimum  field  strengths  recorded 
In  an  hourly  period  for  a  week  are  averaged.  The  average 
difference  over  the  five  week  recording  period  is  then 
determined  from  the  weekly  average.  The  average  difference 
for  the  period  at  the  different  hours  in  shown  for  comparison 
purposes  in  Figure  I5.  The  maximum  average  difference  occurs 
at  the  same  mode  frequencies  as  the  average  values  shewn  in 
Figure  9.  This  s33.:is  to  indicate  that  the  activity  at  the 
mode  frequencies  is  not  only  greater  on  the  average  but  also 


more  varied. 
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Average  difference  of  field  strength  during 
0900-1500  (mst)  for  period  51  July  to 
8  September  I96I 


Conclusions 


Evidence  of  earth* Ionosphere  cavity  mode  frequencies  has 
been  experimentally  obtained.  The  results  Indicate  that  the 
lower  mode  frequencies  are 

fj^ss  8.5  cps 

^2  ~  15.5  cps 
and  f^  ^  22.6  cps. 

The  comparatively  small  amount  of  data  secured  during  the 
evening  hours  precludes  noting  any  change  In  the  mode  frequencies 
due  to  changes  In  Ionosphere  height  or  conductivity. 
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Fort  Trximibull 
New  London ,  Conn . 

Attn:  Mr.  G.  M.  Milligan 

Com.  Officer  and  Director 
US  Naval  Ordnance  Lab. 

Corona.  Calif. 

Attn:  Mr.  A.  W.  Walters 
Code  43 

Technical  Director 

US  Naval  Ordnance  Test  Station 

China  Lake.  California 

Com.  Officer  and  Director 
US  Naval  Ordnance  &  Lab. 

White  Oak.  Maryland 

Attn:  Mr.  Robert  Miller 
Code  042 

Applied  Physics  Laboratory 
Johns  Hopkins  University 
Attn:  Cdr.  Follow 
8621  Georgia  Avenue 
Silver  Springs.  Md. 

Director 

Woods  Hole  Oceanographic  Institution 
Woods  Hole.  Mass.  (Attn.  Dr.  Hersey) 

Dr.  Kenneth  A.  Norton,  Chief 
Radio  Propagation  Bngr.  Div. 

Central  Radio  Prop.  Lab. 

National  Bureau  of  Standards 
Bou Ide r .  Co lor ado 

National  Research  Council 
Com.  on  Undersea  Warfare 
2101  Constitution  Ave..  N.W. 

Wshington  23 >  D.  C. 

Dr.  J.  R.  Wait.  Consultant 
Central  Radio  Prop.  Lab. 

NBS 

Boulder.  Colorado 
Hoffman  Radio 

R.  A.  Schor.  Unit  Supervisor 
Communication  Section 
P.O.  Box  2471 

Los  Angeles  34.  California 


Chief  of  Naval  Research 
Navy  Department 
Washington  23 <  D.  C. 

Code  466 
Code  427 

Ordnaunce  Research  Lab. 

Pennsylvania  State  College 

P.O.  Box  30 

State  College,  Pa. 

Director, 

Scripps  Institution  of  Oceanography 
La  Jolla,  Calif. 

Director , 

Hudson  Labs. 

Columbia  University 

P.O.  Box  329 

Dobbs  Ferry,  New  York 

Dr.  Cullen  Crain 
Rand  Corporation 
Santa  Monica,  California 

University  of  Chicago 
Laboratories  for  Appl.  Sci. 

Museum  of  Science 
Chicago  37 «  Ill. 

Attn:  Mr.  Van  Zeelind 

Commanding  Officer 
Office  of  Naval  Research 
Branch  Office 
1030  East  Green  Street 
Pasadena  1,  Calif. 

National  Science  Foundation 
Engineering  Program 
Washington  23>  D.  C. 

Director 

Research  Lab.  of  Elec. 

Mass.  Institute  of  Technology 
Cambridge  39 >  Mass. 

Dr.  R.  W.  P.  King 
Cruft  Laboratory 
Harvard  University 
Cambridge  38 >  Mass. 

Dr.  R.  H.  Duncan 
Physical  Science  Lab. 

New  Mexico  State  University 

University  Park 

Las  Cruces,  New  Mexico 


P.  S.  Carter 

Radio  Corporation  of  America 
Rocky  Point,  L.  I.,  New  York 


1 


Dr.  J.  H.  Mulligan,  Jr. 

Chairman,  EE  Department 
New  York  University 

New  York  53.  Mew  York  1 

Mr.  Loren  S.  Bearce 
Code  5423 

Naval  Research  Laboratory 

Washington  25.  O*  C.  1 

The  Mitre  Corporation 
Middlesex  Turnpike 
Bedford,  Massachusetts 

Attn:  Mr.  Karl  Swartzel  1 

Dr.  Ronald  V.  Row 

Sylvania  Elec .  Systems 

Div.  of  Sylvania  Elec.  Products,  Inc. 

100  First  Avenue 

Waltham  54,  Massachusetts  1 

Mr.  E.  L.  Sanderson 

US  Naval  Mine  Detection  Laboratory 

Panama  City,  Florida  1 

Dr.  L.  Katz 

Applied  Physics  Laboratory 
Johns  Hopkins  University 

Silver  Springs,  Maryland  1 

Program  Director 
Advanced  Science  Programs 
National  Science  Foundation 
1951  Constitution  Avenue 

Washington  25.  O.  C.  1 


